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ABSTRACT 

This paper presents an analysis of the concentration of the hydrogen molecule inside the 
ionized region of planetary nebulae. The equations corresponding to the ionization and chemical 
equilibria of H, H + , H — , H2, Hjj , and H;j~ are coupled with the equations of ionization and 
thermal balance for a photoionized atomic gas. Forty different reactions related to the formation 
or the destruction of these species are included. The presence of dust is taken into account, since 
grains act as catalysts for the production of H2 , as well as shield the molecules against the stellar 
ionizing radiation. We analyze the effect of the stellar ionizing continuum, as well as of the gas 
and grain properties on the calculated H2 mass. It is shown that a significant concentration of 
H2 can survive inside the ionized region of planetary nebulae, mostly in the inner region of the 
recombination zone. The total H2 to total hydrogen mass ratio inside the ionized region increases 
with the central star temperature, and, depending on the PN physical conditions, it can be of 
the order of ~ 10 ~ 6 or even higher. The increase of the recombination zone with the stellar 
temperature can account for such correlation. This can explain why the H2 emission is more 
frequently observed in bipolar planetary nebulae (Gatley's rule), since this kind of object has 
typically hotter stars. Applying our results for the planetary nebula NGC 6720, we obtain an H 2 
to hydrogen mass ratio similar to the value obtained from the observed H 2 line emission. 

Subject headings: astrochemistry — ISM: molecules — planetary nebulae: general 



1. INTRODUCTION 

Since the first detection of molecular hydro- 
gen in NGC 7027 (Treffers et al. 1976), this 
molecule has been detected in many planetary 
nebulae (PNe), in both the infrared (Aspin et al. 
1993; Kastner et al. 1996; Hora, Latter, & Deutch 
1999) and the ultraviolet (Bowers et al. 1995; Her- 
ald & Bianchi 2002). 

The analysis and interpretation of the infrared 
emission due to hydrogen molecules in PNe has 
been made in the context of: (a) dense neutral 
clumps inside the ionized region (Beckwith, Pcrs- 
son, & Gatley 1978; Reay, Walton, & Atherton 
1988; Tielens 1993; Schild 1995), (b) photodis- 
sociation regions (PDRs; Tielens 1993; Natta & 
Hollenbach 1998; Vicini et al. 1999), as well as (c) 
shocks between the expanding shell and the pre- 
cursor red giant wind (Natta & Hollenbach 1998). 
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Visual inspection of H 2 and [Nil] line images ob- 
served by several authors (Schild 1995; Guerrero et 
al. 2000; Arias et al. 2001; Bohigas 2001) can not 
exclude, however, that the H2 emission is partly 
(or mostly) produced inside the ionized region. In 
fact, the presence of grains in these regions, the 
coexistence of electrons and hydrogen ions, and 
the corresponding thresholds for the destruction or 
ionization of the different species, suggest that H 2 
molecules can survive inside ionized regions. In- 
terpretation of the far ultraviolet absorption spec- 
tra recently obtained by the Far Ultraviolet Spec- 
troscopic Explorer (FUSE) also suggests the ex- 
istence of a hot component of excited H2, which 
could be inside the ionized gas in PNe (Herald & 
Bianchi 2002). 

Formation on grain surface is the main pro- 
cess for the production of H2 molecules in the 
interstellar medium. This formation route was 
first proposed to explain the amount of H2 found 
in interstellar clouds (Hollenbach & Salpeter 
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1971). Grains are condensed in the atmospheres 
of asymptotic giant branch and red giant stars 
(see Whittet 1992 and Evans 1994) and, since 
these stars correspond to the earlier evolution- 
ary stages of PNe, grains could be present in the 
ionized region if they are not destroyed by the 
stellar radiation. Continuum infrared radiation 
observations and modeling of PNe indicate that 
dust is present inside the ionized region of PNe 
(e.g., Natta & Panagia 1981; Harrington, Monk, 
& Clegg 1988; Lenzuni, Natta & Panagia 1989; 
Stasinska & Szczerba 1999; Van Hoof et al. 2000). 
In particular, the distribution of the infrared emis- 
sion follows the morphology of infrared lines and 
continuum radiation emitted by the ionized gas 
(Barlow 1993; Meixner et al. 1996), also suggest- 
ing the presence of grains in the ionized region. 
A very strong evidence that dust is present in the 
ionized region comes from the models obtained by 
Kingdon, Fcrland, & Feibelman (1995) and King- 
don & Ferland (1997), which show depletion of 
elements that can be present in grains. Some au- 
thors (for example, Stasinska 2002 and references 
therein) argue, however, that the grains are in 
clumps and condensations inside the ionized re- 
gion. If grains are indeed mixed with the ionized 
gas, they can act as catalyst for the formation of 
H 2 molecules. 

The presence of electrons and/or H + inside the 
ionized gas can also produce hydrogen molecule 
nuclei, through radiative attachment and associa- 
tive detachment reactions. On the other hand, 
the threshold for the photodissociation of hydro- 
gen molecules is 14.7 eV, while the photoioniza- 
tion threshold is 15.4 eV. These values are higher 
than the threshold for the ionization of hydrogen 
atoms (13.6 eV). Hydrogen molecules could then 
be shielded from the photons that could dissociate 
it, surviving inside the ionized region in zones far 
from the ionizing source. 

Detailed analyses of the presence of molecules 
have been done for interstellar clouds (for exam- 
ple, de Jong 1972; Black & Dalgarno 1977; Viala, 
Roueff, & Abgrall 1988), as well as for photodisso- 
ciation regions (for example, Ticlens & Hollenbach 
1985; Draine & Bertoldi 1996). Models for pho- 
todissociation regions have been applied to plane- 
tary nebulae, assuming that the incident radiation 
is filtered after crossing the ionized region (e.g., 
Natta & Hollenbach 1998). These models assume 



that the flux of the incidentc radiation is zero be- 
tween 13.6 and 54.4 eV. Due to the behavior of the 
ionizing cross-sections with the photon energy, the 
correct ionizing spectrum that reaches each point 
of the ionized region is, however, much more com- 
plex, and a proper account of the radiation trans- 
fer along the nebula must be done when analyzing 
the molecular chemistry. 

Since the first detection of molecular hydrogen 
in PNe, only a few attempts to calculate the H 2 
content in ionized regions have been made (Black 
1978; Gussie & Pritchet 1988; Cecchi-Pestellini & 
Dalgarno 1993). In these studies, however, one or 
more of the following effects is not taken into ac- 
count: the effect of dust on the radiation transfer, 
the formation of H 2 molecules on grain surface, 
and the photodissociation in two steps. 

An interesting correlation has been obtained 
between the presence of H 2 line emission and the 
morphological type of the planetary nebula. Sig- 
nificant H 2 emission is more frequently observed 
in PNe having bipolar morphology (Gatley's rule; 
Zuckerman & Gatley 1988; Kastner at al. 1996). 
This correlation between H 2 emission and the mor- 
phology is not understood yet. 

Furthermore, some problems regarding the line 
emission of H 2 and neutral atoms in PNe are not 
solved yet. The excitation mechanism of the ob- 
served H 2 emission lines, for example, is still un- 
certain. Some authors conclude that excitation is 
due to radiative cascades following the UV pump- 
ing. Shock waves produced by the expanding 
nebulae are also suggested as a possible excita- 
tion mechanism. See, for example, Black, Porter 
& Dalgarno (1981), Schild (1995), Kastner et al. 
(1996), Natta & Hollenbach (1998), Hora et al. 
(1999), Guerrero et al. (2000), and Arias et al. 
(2001). Another problem is that the emission 
lines produced in low-ionization regions by neu- 
tral atoms like C, N, and O are not well fitted 
by the models (Aller 1984; Clegg et al. 1987; van 
Hoof et al. 2000; Stasinska & Szczerba 2001). The 
analysis of these lines is made using photoioniza- 
tion models, where lines of neutral atoms are pro- 
duced in the outer regions of a photoionized gas, 
and the presence of molecules is not taken into 
account. The hydrogen molecule, being the most 
abundant, can affect the physical conditions of the 
gas, since its presence changes the concentration 
of the different species in the gas, as well as the 
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thermal equilibrium, due to the mechanisms in- 
volving H 2 . If molecules are indeed present in the 
region where lines of neutral C, N, and O are pro- 
duced, and if their effect on the gas is not taken 
into account, the results can be misleading. A de- 
tailed study, with a self-consistent modeling of the 
presence of H2 molecules, is then important to a 
better knowledge of the chemical composition and 
physical conditions of PNc. 

The aim of this paper is to make a first estimate 
of the amount of the H2 molecule in the ionized 
regions of planetary nebulae, in order to verify if 
the presence of this molecule can be significant in 
these regions, as suggested by the observations. If 
the calculations show that the amount of H 2 is im- 
portant, a next step would be a careful calculation 
of the level population distribution, in order to 
properly account for the different mechanisms that 
depend on the level population. In the present pa- 
per, the equations for the chemical equilibrium of 
molecular hydrogen are coupled to the equations 
of a photoionization code. The concentration of 
the hydrogen molecule is obtained along the neb- 
ula, as well as the total amount of molecules inside 
it. The theoretical models used in our analysis are 
described in § 2. The results for a standard PN 
and the effect of the input parameters on the ob- 
tained H 2 mass appear in § 3, while a discussion 
of the main results is presented in § 4. 

2. MODELS 

Models for typical PNe are obtained with the 
photoionization code AANGABA (Gruenwald & 
Viegas 1992). The physical conditions of the gas 
are determined by solving the coupled equations 
of ionization and thermal balance for a spheri- 
cally symmetric cloud. Several processes of ion- 
ization and recombination, as well as of gas heat- 
ing and cooling, are taken into account. The en- 
ergy input to the gas is essentially due to ioniza- 
tion of atoms and atomic ions by primary (stel- 
lar) and secondary (diffuse) radiation. The pro- 
cesses contributing to the gas cooling are colisional 
excitation, radiative and diclectronic recombina- 
tion, thermal ionization and free-free emission. 
The transfer of the primary and diffuse radiation 
fields is treated in the outward-only approxima- 
tion. Twelve elements (H, He, C, N, O, Mg, Ne, 
Si, S, Ar, CI, and Fe) are included. The code was 



modified in order to include the chemical equilib- 
rium of the species H, H+, H , H 2 , , and H3 . 
The corresponding equations were then coupled 
with the equations for the ionization balance for 
the atomic gas. The H 3 molecule is not included 
due to its short lifetime (~ 10~ 7 s; Gaillard et al. 
1983; Bordas, Cosby, & Helm 1990). The ther- 
mal balance was solved for the atoms and atomic 
ions; the effect due to H~, H 2 , rLj", and H3 on the 
thermal balance is not taken into account. 

Several processes that can form or destroy the 
included species are taken into account, as for ex- 
ample chemical reactions, photoprocesses, as well 
as grain surface reaction. The forty reactions in- 
cluded in the hydrogen equilibrium are listed in 
Table 1. Rate coefficients and cross sections are 
taken from the literature (sec references in Table 
1). Special attention was given to the search of the 
reactions and their corresponding rate coefficients, 
in order to have a set as complete and accurate as 
possible for the physical conditions prevailing in 
PNe. Some rate coefficients depend on the distri- 
bution of the reactants on their different energy 
levels. As said above, we do not calculate in de- 
tail the level distribution of the molecules in the 
present paper; a LTE level distribution for the hy- 
drogen molecules, taking into account the three 
lowest vibrational levels of the electronic ground 
state is, however, adopted when cross-sections for 
the photoprocess are provided by the literature for 
individual levels. 

Grains are assumed to be uniformly mixed with 
the gas, and are heated by the stellar and nebu- 
lar continuum radiation. The attenuation of the 
stellar and diffuse radiation by dust is taken into 
account; the effect of scattering by dust on the 
radiation field is, however, not included. 

2.1. The H 2 Formation on Grain Surface 

The formation of H 2 on grain surface requires 
an analysis of the grain survival. We assume that 
the grain can survive if the grain temperature 
(T g ) is smaller than the sublimation temperature, 
which is a characteristic of the grain material. 

The grain temperature is calculated assuming 
radiative equilibrium. The grains are heated by 
the stellar and the nebular continuum radiation, 
and emit as blackbodies. 

Following Ticlcns & Hollcnbach (1985) and Hol- 
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lenbach & Salpeter (1971), the rate of H2 forma- 
tion on grain surface, T g , is given by the product of 
the rate of H-grain collisions by probabilistic fac- 
tors that depend on the grain surface properties 
and on the gas temperature, i.e., 

S f 

T g = —%a 2 g v H n g n{H) (1) 

where a g is the grain radius, vr is the thermal ve- 
locity of the H atoms, n g is the grain density, and 
n(H) is the neutral atomic H density. The prob- 
abilistic factors S and / are given by Hollenbach 
& McKee (1979). These factors decrease with the 
grain temperature (T g ) and the gas temperature 
(T), particularly when T g is higher than 100K. 
They also depend on the properties of the grain 
material. 

A new model for the H2 formation in grains is 
suggested by Cazaux & Tielens (2002). In this 
model, they include the physiobsorption of the 
H atom (which dominates in temperatures lower 
than 50-100K), as well as the chemisorption, which 
would be more important for higher temperatures. 
However, as stated by Glover (2003), the subject 
is still not clear. 

From equation (1), the rate of H2 formation 
on grain surface is proportional to the total avail- 
able surface area of all grains, i.e., Y g oc a g n g . 
However, following the most common notation in 
the literature, we chose to describe the quantity of 
grains by means of the dust-to-gas ratio, 

M d = ^a 3 g p g n g 

M g 3m H n H (l + ^n H e/n H ) ' 

where tuh is the atomic mass of H, nn and 
are the total numerical density of hydrogen and 
helium nuclei, respectively, and p g is the grain spe- 
cific density. Assuming a Maxwellian distribution 
of velocity for the H atoms, equation (1) can then 
be rewritten as 

SfT-^M d /M g n H n{H) 

We define for convenience (M d /M 9 )a~ 1 hereafter 
as [3. 

2.2. Model Input Parameters 

A grid of theoretical models is obtained with 
input parameters typical of PNe, which includes 



the incident spectrum, the gas density and the el- 
emental abundances of the nebula, as well as the 
dust material, size, and density. 

In our models, a blackbody radiation is as- 
sumed for the PN central star, characterized by 
its temperature (T*) and luminosity (£*). Accord- 
ing to Pottasch (1984), T* ranges from 3xl0 4 to 
5xl0 5 K, while L* ranges from 1O 2 L to 1O 4 L . 

The nebular gas is characterized by its total 
hydrogen nuclei density (n#), which is assumed 
constant inside the nebula, and by its elemental 
abundances. Gas densities typically range from 
10 2 to 10 5 cm~ 3 (Pottasch, 1984). We assume 
constant elemental abundances for the atomic gas 
along the nebulae, with the same values for all 
models. These are mean values for planetary neb- 
ulae taken from Kingsbourgh & Barlow (1994), for 
He, C, N, O, Ne, S, and Ar, obtained from atomic 
line-emission intensities. For the elements whose 
abundance is not provided by these authors (Mg, 
Si, CI, and Fe), we adopted the values given by 
Stasinska & Tylenda (1986), in order to make a 
rough correction for grain depletion. Since the ef- 
fect of these elements is not very effective for the 
gas cooling, their exact proportion in the form of 
grains will not be important. 

Typical grain radius ranges from 10~ 3 /zm to 10 
yitm and dust-to-gas mass ratio from 10~ 5 to 10 _1 
in PNe (Stasinska & Szczerba 1999). The compo- 
sition of the ISM grains (Allcr 1984; Evans 1994), 
i.e., silicate and carbon-based compounds, is as- 
sumed here for the PNe grains. In this work, 
we analyze the presence of amorphous carbon, 
graphite, silicate, smoothed ultraviolet silicate, sil- 
icon carbide (SiC), as well as neutral and ionized 
polycyclic aromatic hydrocarbons (PAHs). Op- 
tical properties for these astrophysical grains are 
taken from Rowan-Robinson (1986), Draine & Lee 
(1984), Laor & Draine (1993), Weingartner & 
Draine (2001), and Li & Draine (2001). These 
properties depend on the grain size and composi- 
tion, as well as on the photon energy. Sublimation 
temperatures and specific densities are taken from 
Laor & Draine (1993) and Li & Draine (2000). 

We define a standard PN with a given set of 
input parameters, in order to explore the influ- 
ence due to each one on the produced amount of 
H2. For this, we vary each parameter along its 
typical range, while keeping the other parameters 
fixed. The input parameters for the standard PN 
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are the following: T* = 150000 K, L* = 3OOOL , 
n# = 10 3 cur 3 , a 3 = 0.1/zm, M d /M g = 10~ 3 , 
and amorphous carbon grains. 

3. RESULTS 

3.1. Standard Model 

The calculated concentration of species contain- 
ing hydrogen nuclei versus the distance from the 
central star is shown on the upper panel of Figure 
1, for the standard PN model. The concentrations 
are given in units of uh ■ The temperature radial 
profile for the standard nebula is plotted on the 
lower panel. The distance is given in units of the 
outer radius of the ionized region, R ma x, which is 
defined as the distance to the central star where 
the ionization degree, n(H + )/nn, reaches 10 -4 . 
For the standard PN, R ma x — 0.21 pc; for a neb- 
ula with same input parameters, but T* = 50000 
K, R m ax = 0.17 pc, and, if the nebula has a den- 
sity tih = 10 4 cm~ 3 and other parameters as for 
the standard nebula, R ma x = 0.04 pc. 

As it can be seen in Figure 1, the hydrogen 
molecule is mainly ionized in the inner and in- 
termediate regions. The abundance of the neutral 
molecule, as well as of the ionized molecule, has its 
maximum in the inner region of the recombination 
zone, decreasing slightly in the less ionized regions. 
This behavior can be explained by the curves of 
Figure 2, which shows the rates for the most im- 
portant reactions of formation and destruction of 
the hydrogen molecule in the standard PN, as a 
function of the distance from the central star. 

In brief, the most important reactions for the 
formation of the hydrogen molecule nuclei (H2 or 
H2 ) are: 

H + H~ -> H 2 + e~ (4a) 
2H + grain — > H 2 + grain (4b) 
H + H+ ->H+ + hv (4c) 
H+ + H- -► H+ + e~ (4d) 

with H~ produced by the reaction 

H + e~ H~ + hv (5a) 

Charge exchange reactions are effective in ionizing 
or recombining the hydrogen molecule: 

H+ + H^H 2 + H+ (6a) 



H 2 +H+ -> H+ + H (6b) 

(6c) 

Comparing the distribution of the various 
species (Fig. 1) with the formation and destruc- 
tion rates along the standard PN (Fig. 2), the 
maximum concentration of molecular nuclei can 
be associated to a semi-ionized region, where H, 
H + and c~ coexist with significant abundance. In 
fact, as said above, the maximum abundance of 
molecular hydrogen occurs in the inner recombi- 
nation zone, and the distribution follows a radial 
distribution similar to that of H . 

Grains are the most important route to the for- 
mation of H 2 in the outer regions since other re- 
actions include ionized species. Therefore, when 
grains are not present the concentration of the 
molecule decreases in the outer regions. 

Regarding the destruction mechanisms, pho- 
todissociation is an important process for the de- 
struction of H 2 , as well as of H^. In particular, 
the photodissociation in two steps is important for 
the destruction of H 2 in the whole nebula, since 
the maximum threshold for this process is lower 
than that of hydrogen. The direct photodissocia- 
tion, as well as photoionization of H2 are of minor 
importance, mainly in the intermediate and outer 
zones, since the energy threshold is higher than 
13.6 eV. 

The abundance of is much lower compared 
with that of other species, and has its maximum in 
the same region as H 2 and H 2 ■ The formation of 
this molecular ion occurs mainly through radiative 
association of H 2 with H+, but the ion- molecule 
reaction of H 2 and R 2 ^ s & ^ so important in the re- 
combination zone. Dissociative recombination and 
collisions with H are the main destruction routes 
ofH+ 

3.2. The Effect of Dust 

The grain temperature depends on the incident 
radiation, as well as on the grain optical proper- 
ties, which, in turn, depend on the grain material 
and size. The grain temperature decreases steeply 
in the inner 20% region of the nebula, and de- 
creases slowly for increasing radius. For example, 
assuming a g < 0.0 1/mi, and all other input param- 
eters the same as those for the standard model, the 
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temperature for graphite grains decreases from ~ 
500 K, in the inner regions, to ~ 50 K, in the far- 
thest zones, while for a g < 0.1/zm, the grain tem- 
perature ranges from 180 K to 26 K. For other ma- 
terials the temperatures are slightly lower. Grains 
can survive inside the entire ionized region, except 
those made of PAH, which sublimate for distances 
shorter than <~ 0.01R max . 

Direct stellar radiation is much more important 
for the grain heating than absorption of Lya radi- 
ation in the ionized region. The grain temperature 
is increased by ~ 10% when Lya radiation is in- 
cluded in the grain heating. Such increase leads 
to a maximum decrease of only 0.3% in the H 2 
concentration. 

The effect of the dust-to-gas ratio on the rela- 
tive abundance of H2 molecules is shown in Fig- 
ure 3. The results presented in this figure corre- 
spond to models with the same input parameters 
as in the standard model, except for the dust-to- 
gas ratio. In the figure, the behavior of n(7?2) with 
R/Rmax is shown for different values of Md/M g , 
with the lowest curve corresponding to a model 
without dust. Since the reaction on grain surface 
is an important reaction of H2 formation all along 
the nebula, increasing Md/M g results in higher 
values of n(_ff 2 It can be also noticed that 

for R/R m ax between 0.8 and 0.9 the results do 
not depend on Md/M g . This can be explained by 
the fact that in this region the dominant processes 
of H2 formation are the charge exchange and the 
associative detachment reactions. As said above, 
when few or no grains are present, there is a steep 
decrease of H2 molecules in the outer regions of the 
nebula, due to the decrease of H~ , which is mainly 
produced in regions where the electron density is 
significant. 

As discussed in §2.1, the key parameter for the 
formation of H 2 in the surface of grains is the to- 
tal available surface area of grains or, in terms of 
our input parameters, the dust-to-gas mass ratio 
divided by the grain radius (/3) . For a given dust- 
to-gas ratio, smaller grains imply a higher forma- 
tion rate (see eq.[ 3]) of H2 molecules. On the 
other hand, smaller grains have larger tempera- 
tures, leading to lower formation rates. However, 
the effect due to the grain temperature is small; 
fixing the dust-to-gas ratio, a decrease of the size 
of the grain leads to an increase of the abundance 
of H 2 . 



As said in §1, the formation of H2 molecules on 
the surface of grains has been assumed in the in- 
terstellar medium in order to explain the required 
amount of H2 molecules. For the ionized region 
of PNe, however, the associative detachment re- 
action, as well as the charge transfer of H^ with 
H are also important processes for the formation 
of H 2 molecules. The dust contributes more sig- 
nificantly for the production of H2 molecules for 
values of (3 > 100 cm -1 . 

3.3. Other models 

The distribution of the various species, as well 
as that of the reactions rates, along the nebula for 
models with input parameters other than those 
of the standard model is similar to that shown 
in Figures 1 and 2. In all models, the maximum 
H 2 concentration occurs where the gas tempera- 
ture is around 7500 K and the ionization degree 
n(H + ) /nji is approximately 0.1. 

In order to compare the amount of H2 produced 
in different PNe conditions, we define the param- 
eter Rm, which corresponds to the ratio between 
the total hydrogen molecule nuclei mass (H 2 + 
Hj ) and the total hydrogen nuclei mass, i.e., 



R 



M 



J*— An[n(H 2 ) 



n(H2)]mu 2 r 2 dr 



J R max AirnH'mHr 2 dr 



(7) 



where r is the distance from the central star, while 
Rq and Rmax are the inner and the outer radius of 
the ionized region as defined above, respectively. 
The inner radius is chosen as i?o = 10 15 cm, but 
its exact value does not affect the results presented 
here, provided that it is lower than §.SR max . 

For the standard model, we obtain a value for 
Rm of the order of 2xl0~ 6 . Increasing the dust- 
to-gas ratio (or decreasing the grain size) by a fac- 
tor of 10, and keeping the other input parameters 
fixed, the fraction of molecular mass increases by 
a factor around 2. As discussed above, the main 
parameter regarding the presence of grains is the 
available total surface of grains. When dust is not 
present in the gas, Rm is only 10% lower than in 
the standard model. In fact, Rm increases with 
Md/M g and/or decreasing a gi but, as said above, 
for a value of (3 lower than that for the standard 
model, the presence of grains does not affect the 
obtained results for Rm- 
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Different grain materials leads to a maximum 
difference in Rm by less than 25%. 

Plots of Rm versus T» are shown in Figure 
4, for different values of L», uh, and Md/M g . 
In this figure, it can be seen that Rm increases 
with T*. This is an interesting result, since hot- 
ter stars produce more high-energy photons. How- 
ever, since the hydrogen photoionization cross sec- 
tion decreases with the photon energy, high-energy 
photons can travel farther into the nebula before 
being absorbed by the gas, producing a more ex- 
tended recombination zone. The increase in Rm 
with T» can explain why H2 is more frequently 
detected in bipolar PNe, since these objects have 
hotter central stars (Corradi & Schwarz 1995). In 
fact, the relative thickness of the recombination 
zone is the crucial factor for Rm, and it implies 
that Rm also increases with the decrease of n# 
and Z/*, as shown in Figure 4. An increase of the 
gas density by a factor 10 decreases Rm by a fac- 
tor of the order of 3, while a change of a factor 
100 in the luminosity leads to a change of Rm by 
a factor around 24. However, Rm increases more 
steeply with T* than with any other input param- 
eter, as can be seen from Figure 4. An increase of 
T* by a factor 10 leads to an increase of Rm by a 
factor around 6000. 

4. CONCLUDING REMARKS 

Using a photoionization code which includes 
atoms and ions of twelve elements, as well as H , 
H 2 , H^, and H3 , we analyze the presence of the H 2 
molecule over the entire ionized region of PNe, up 
to the region where the ionization degree decreases 
to 10 -4 . The transfer of the radiation emitted by 
the star, as well as of that produced by the neb- 
ula, is treated in detail in order to obtain a self- 
consistent distribution of the hydrogen molecule in 
the nebula. The presence of grains is taken into ac- 
count, since they are present in the ionized region 
and can act as catalyst in the formation of the H 2 
molecule. The effect of the various reactions that 
contribute to the formation and destruction of H 2 
is studied. 

Our calculations show that H 2 molecules can 
survive inside the ionized region of PNe, particu- 
larly in the inner zone of the recombination region, 
where the ionization degree is of the order of 0.1. 
In this region, the formation rate of H 2 molecules 



is enhanced due to the high efficiency of the as- 
sociative detachment of H and H~ and the charge 
exchange between H and H^, while the destruction 
is not very effective due to the shielding provided 
by H atoms. When significant superficial area of 
grains is available (i.e., j3 > 100 cm -1 ), the grain 
reaction is the main process of H 2 production in 
regions of low ionization. 

Notice that the grain surface reaction is not 
the only important production mechanism of H 2 
molecules in PNe, as it is in the interstellar 
medium. The associative detachment of H and 
H~ , as well as the charge exchange between H and 
Hj are also important processes for the formation 
of H 2 molecules. In the H recombination zone, 
where the density of H 2 reaches its higher values, 
the main processes, even when large amounts of 
grains are present, are via associative detachment 
and charge exchange. The H and the H^" ions 
have similar concentrations, except in the outer 
region of the PNe, and are abundant in regions 
where both electrons and H ions are present. 

The effect of the stellar temperature and lumi- 
nosity, as well as that of the gas density, dust- 
to-gas mass fraction, and grain material, on the 
obtained abundance of the hydrogen molecule, is 
also analyzed. Regarding the presence of grains, 
the available superficial area is the main parame- 
ter that controls the molecule formation efficiency, 
while the grain temperature and composition are 
of minor importance. 

We conclude that Rm, the H 2 to total hydrogen 
nuclei mass ratio, increases with the relative thick- 
ness of the recombination zone. This explains why 
Rm increases with T* and decreases with nn and 
£*, the main effect being due to the stellar temper- 
ature. The behavior of Rm versus T* can explain 
why the emission of H 2 is more often detected in 
PNe with bipolar morphology, since these objects 
have in general hotter central stars. 

The presence of molecular material is gener- 
ally used as evidence that the region is optically 
thick. However, if the hydrogen molecule can sur- 
vive within the hot gas, the detection of H 2 emis- 
sion is not an indication that a nebula is optically 
thick. This point has important consequences for 
the determination of PNe distances, since some 
distance indicators are obtained assuming that the 
gas is completely optically thick (i.e., all ionizing 
photons are absorbed). The same applies to the 
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Zanstra method for the determination of the stel- 
lar temperature. 

As shown in this paper, H 2 molecules can sur- 
vive not just in very low, but also in mild ioniza- 
tion regions inside PNe. Thus, a detailed study of 
the effect the presence of the hydrogen molecule 
on the gas thermal equilibrium is essential for an 
accurate picture of the physical processes inside 
PNe. The presence of H2 molecules in mild ionized 
gas can solve, for example, some problems faced 
by photoionization models in explaining emission 
lines from the low-ionization regions. Hydrogen 
molecules can affect the physical conditions of the 
gas, providing the additional heating required to 
explain the discrepancies. A more realistic esti- 
mate for the gas abundance could also be obtained. 

In this work we do not include the self-shielding 
of H2 in the calculations, since it would be nec- 
essary to obtain the distribution of the molecule 
in its various energy levels in order to estimate 
the shielding. The effect of the self-shielding on 
the photoprocessess can be very important in the 
outer regions of PNe where the H2 column density 
is high. Approximated methods for the treatment 
of the self-shielding effect on the photodissocia- 
tion of H2 are available in the literature as, for 
example, in Van der Werf (1989) and Draine & 
Bertoldi (1996). Taking this approximated meth- 
ods into account, the computed H2 concentration 
can increase up to two orders of magnitude in the 
outer parts of the ionized region. Therefore the 
calculated H2 concentrations in the present work 
can be considered as lower limits. 

Shocks arc frequently invoked in the literature 
for explaining the H2 line emission since pure ra- 
diative excitation does not account for the obser- 
vations. The presence of H 2 molecules in regions 
of higher excitation level, excited by UV radiation 
not fully absorbed by inner regions, could account 
for the line emission, without invoking shocks. 

Molecular and atomic hydrogen masses for the 
Ring Nebula (NGC 6720) were derived from ob- 
servations by Greenhouse, Hay ward, & Thronson 
(1988), resulting in a molecular to atomic gas ratio 
equal to ~ 7 x 1CT 5 . Assuming the equilibrium be- 
tween formation on grain surface and photodisso- 
ciation, these authors conclude that chemical equi- 
librium does not account for the detected quantity 
of H 2 . However, from our more detailed model and 
assuming the parameters for the central star and 



the gas given by Greenhouse et al. (1988), we ob- 
tain Rm = 2.4 x 10 -5 , which is reasonably similar 
to the value obtained by these authors. 
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tions for the hydrogen molecule, and to T. Abel 
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are pleased to thank the anonymous referee for his 
valuable comments and suggestions that helped to 
improve the paper. The authors thank FAPESP, 
CNPq, and Pronex/CNPq for the financial sup- 
port of this work. 
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R/R, 



Fig. 1 . — Concentrations of the species containing 
H nuclei (upper panel) and gas temperature (lower 
panel), for the standard model. 



0.0 0.2 0.4 0.6 0.8 1.0 X 0.0 0.2 0.4 0.6 0.8 1.0 
R /R m0 . R / R ™» 

Fig. 2. — Rates of the most important formation 
and destruction reactions of molecular hydrogen, 
for the standard model. H2 formation (upper left): 
H + H~ (solid curve), + H (dashed), and H 
+ H + grain (dot-dashed). H 2 destruction (upper 
right): H 2 + hv (photodissociation in two steps; 
solid), H2 + hv (direct photodissociation; dashed), 
H2 + H + (charge exchange, dot-dashed) , H2 + e~ 
(colisional dissociation, dot-dot-dot-dashed), and 
H 2 + hv (photoionization; dotted). formation 
(lower left): H 2 + H+ (solid), H + H+ (dashed), 
and H + H + (dot-dashed). destruction (lower 
right): + hv (producing H + + H; solid), 
+ e~ (dashed), + H (dot-dashed), and + 
hv (producing 2H + + e~; dotted). 
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Fig. 3. — Effect of the dust-to-gas mass ratio on 
the calculated abundance of H2. Curves are shown 
for M d /M g equal to 0, 1CT 3 , 1(T 2 , and 10" 1 ; 
the lower curve corresponds to a model without 
grains, while upper curves correspond to models 
with higher values of the dust-to-gas ratio. Other 
parameters are the same as in the standard PN 
model. 
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Fig. 4. — Molecular hydrogen mass fraction ver- 
sus central stellar temperature. Upper panel: The 
results are shown for nn = 10 3 , 10 4 , and 10 5 
cm -3 (lower curves correspond to higher nn) for 
two values of (3: 10 2 cm -1 (solid curves) and 10 4 
cm -1 (dashed). Lower panel: Results for L* = 
500, 3000, and 10 4 L Q , with Rm decreasing with 
increasing L*. Other parameters correspond to 
those of the standard model. 
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Table 1 

Chemical Reactions and Physical Processes 



Reaction 



Ref. 



Reaction 



Ref. 



Photoionization 
H + hp -> H+ + e" 
H 2 + hv -> + e~ 

Photodetachment 
H~ + hv ^ H + e~ 

Radiative Recombination 
ff + + e- -> i? + 

Dissociative Recombination 



2H 
H 2 - 
3H 



H 



Photodissociation 
H 2 + hv -> 2ff 
H 2 + hv 

H + H+ 
2H+ + e 



H+ + hv 
H+ + hv 



Radiative Attachment 
H + er ^ H- + hv 



Charge Exchange 
H+ + H ^H+ + H 2 
H 2 + H+ + H+ 
He+ + H -> He + H+ 
He++ + H -> He+ + H+ 



12 
13 
14 
15 



16 



17 
7 
18 
19 



Associative Detachment 
H + H- ^H 2 + e~ 
H++H- -> JT+ + e" 



Radiative Association 



H + H+ - 
H 2 + H+ 



hv 



H+ + hv 



Ion-Molecule Reaction 
H 2 + +H 2 ^H+ + H 



Grain Surface Reaction 
2if + grain — > i? 2 + grain 



20 



Neutralization 



Other Processes 



H+ +H- 


-» 2i? 


9 


H 2 + e" -. 


■ 2H + e~ 


21 


H+ + H~ 




9, 10 


ff 2 + e- - 


■ H + H~ 


22 


H+ + H~ 


-» 3i? 


9, 10 


H 2 + H^ 


3H 


23 


H+ + H- 


^2ff 2 


9, 10 


2H 2 — > i?2 


+ 2H 


23 


H+ + H- 


-» i? 2 + 2i7 


9, 10 


ff- +e- - 


^H + 2e- 


16 








H- +H - 


>2H + e- 


22 


Colisional Ionization 




H+ + H 2 - 


-^H + H+ + H 2 


11 


H + er ^ 


H+ + 2e- 


9 


H+ + H- 


> H 2 + H+ 


7 


2H -> H + H+ + e~ 


11 


H+ + H 2 - 


-+H + H 2 + H+ 


11 


2H ^ H + H* ^ H + H+ + e- 


11 


H+ + H 2 - 


-^H+ + 2H 2 


11 
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(5) Schneider et al. (1994, 1997); (6) McCall et al. (2003), with branching 
ratios given by Datz et al. (1995); (7) Galli & Palla (1998); (8) Thcard & 
Huntress (1974); (9) Lepp, Standi, & Dalgarno (2002); (10) following Dalgarno 
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2H; (11) Hollcnbach & McKee (1989); (12) Allison & Dalgarno (1969), with LTE 
distribution for the H2 energy levels; (13) Abgrall et al. (1994) and Abgrall, 
Roueff, & Drira (2000), with LTE distribution for the H 2 energy levels; (14) 
Stancil (1994), with LTE fitting to T = 8400 K; (15) Bates & Opick (1968); 
(16) Abel et al. (1997); (17) Karpas, Anicich, & Huntress (1979); (18) Jura & 
Dalgarno (1971); (19) Arthurs & Hyslop (1957); (20) see text; (21) Donahue & 
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